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COMPARISON OF MODELED AND MEASURED DEPOSITION 
FOR THE WITHLACOOCHEE SPRAY TRIALS? 


J. E. Rafferty”, R. K. Dumbauld? 
H. W. Flake, dr.*, J. Ww. Barry* and J. Wong* 


ABSTRACT 


The Wrthlacoochee spray trials were planned and conducted as a ptlot 
project by Southeastern Area-Forest Pest Management (SA-FPM, and FPM-Methods 
Application Group to evaluate the feastbiltty of ustng aertal spray appltca- 
tton techniques tn the control of cone losses due to coneworm infestations 
tn ptne seed orchards. As part of the overall ptlot project, deposttton 
measurements were made at the orchard canopy top, within sample trees and 
on the ground beneath the trees. These measurements were made for use tn 
defining the effectiveness of the aertal spray appltcatton technique and to 
evaluate the potenttal usefulness of the Cramer-Barry-Grim (CBG) computertzed 
model tn the planning, conduct and analysts of seed orchard spray projects. 
Thts report descrtbes the Wtthlacoochee spray trtals and compares the 
deposttton measurements with the deposttton estimates calculated using the 
CBG model computer program. 


INTRODUCTION 


1. Background 


The USDA-Forest Service (FS) has been using mathematical models to 
predict the atmospheric transport, dispersion and deposition of aerial ly- 
applied pesticides for more than a decade. Simplified aerial line source 
models developed for the US Army were applied (GCA Corporation, 1971; 

H. E. Cramer Company, 1973) early in the decade to determine optimum swath 
widths and application rates for use in pilot tests of insecticides under 
consideration at that time for control applications in western forests. 
The implications of these early efforts in the use of mathematical models 
to improve the planning, conduct and analysis of spray programs were re- 
viewed in a paper presented at the USFS Workshop for Aerial Application 


Prepared for USDA Forest Service, Forest Pest Management, Methods 
Application Group, Davis, CA under Contract No. 53-91S8-9-6260; 
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H. E. Cramer Company, Inc., P. 0. Box 8049, Salt Lake City, UT 84108. 
USDA Forest Service, SA-FPM, Asheville, NC 
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of Insecticides Against Forest Defoliators, held in Missoula, Montana, 
23-29 April 1974 (Dumbauld, Cramer and Barry, 1975 pReatner he ©. Cramer 
Company (Dumbauld and Bjorklund, 1977) used aerial spray dispersion models 
to assist the State of Maine, Bureau of Forestry, in determining offset 
distances required for various aircraft to ensure that drift from spray 
blocks in the Maine 1977 spray program posed no environmental hazard to 
exclusion areas (waterways, homes, etc.) in the vicinity of the spray 
blocks. Under the sponsorship of the USDA Expanded Douglas-fir- Tussock 
Moth Research and Development Program, work began in early 1977 on the 
development of a technical data base for use in applying mathematical dis- 
persion models to FS spray projects and the refinement and adaptation of 
existing models to predict spray behavior above and within forest cano- 
pies. This work (Dumbauld, Rafferty and Bjorklund, 1977) resulted in the 
development of the CBG computerized spray dispersion model and the com- 
parison of model predictions with measurements made during selected FS 
spray programs. The CBG computer based model jis comprised of two major 
parts. The first part simulates the deposition of spray drops at the top 
of a forest canopy. The second part, which is based on a model developed 
by Grim and Barry (1975), simulates the penetration of spray drops into 
the forest canopy. 


FPM-MAG decided to contract with the H. E. Cramer Company to conduct 
a joint field effort to demonstrate and advance FS modeling capabilities 
in predicting requirements for delivering effective doses to target pests 
and control spray drift beyond target areas. Shortly after work started 
under the contract, the FPM office at Asheville, NC, proposed a feasibility 
study of aerial application techniques for the control of seed losses due 
to seed and cone insects in pine seed orchards. FPM-Asheville and FPM-MAG 
jointly decided to conduct a pilot project to evaluate aerial application 
technology for southern seed orchards and to use this pilot project as a 
field demonstration for the CBG modeling capability. -The Withlacoochee 
State Seed Orchard located near Brooksville, Florida, was selected in 
December 1979 as the site of the pilot project which was scheduled to be 
conducted during February 1980. As part of the field demonstration pro- 
ject under Contract No. 53-91S8-9-6260, the H. E. Cramer Company, Inc. 
used the CBG model to calculate optimum swath widths and release heights 
for two aircraft types (fixed-wing and helicopter) with spray releases 
during the early- and late-morning hours at the Withlacoochee State Seed 
Orchard (Rafferty and Dumbauld, 1980). H. E. Cramer Company personnel 
also participated in the conduct of the aerial spray trials and were as- 
signed the primary task of assisting the FS in making meteorological 
measurements required for input to the CBG model. 


Aerial pesticide applications to orchards are usually accomplished 
with aircraft flying immediately above the canopy and calibrated to 
Spray at a rate of from 5 to 10 gallons per acre. Thus, spray strategy 
differs markedly from that used in spraying western forests where aircraft 
fly 50 feet above the canopy and apply spray at rates of one gallon or less 
per acre. The CBG model was developed primarily for use with the latter 


strategy and does not incorporate an adequate wake model for use with low- 
flying aircraft. For this reason close agreement between the model and ob- 
served spray deposition pattern was not expected, but the opportunity to 
develop a data base for future application of the CBG model to seed orchard 
Spray activities was welcomed. 


foc Punpose 


The purpose of this report is to describe the final phase of the field 
demonstration of the CBG model. The major topics discussed in the report 
are a comparison of model deposition estimates with deposition measurements 
made during the Withlacoochee spray trials and an overall evaluation of the 
usefulness of the model in planning and conducting orchard spray projects. 


1.3 Report Organization 


Section 2 contains a general description of the Withlacoochee spray 
trials. The inputs required by the CBG model for dispersion, deposition 
and canopy penetration calculations and the values assigned to these inputs 
are described in Section 3. Section 4 contains a description of the pro- 
cedures used in the model calculations and a comparison of the results of 
the model calculations with the measurements made during the spray trials. 
The results are summarized in Section 5 which also contains recommendations 
for further study of the trial data and. for model improvements. A brief 
description of the model used in this study is provided in Appendix A. 
Appendix B contains drop-size distribution data for the selected trials 
used in the coOmparison of model deposition estimates with deposition measure- 
ments. Appendix C presents graphs of observed and model deposition patterns 
for selected trials. 
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SECTION 2 
WITHLACOOCHEE SPRAY TRIALS 


2.1 General 


The Withlacoochee State Seed Orchard is located a few miles north 
and west from the intersection of Highway 98 and Interstate 75 near 
Brooksville, Florida. The surrounding terrain is gently rolling with 
elevations varying between 50 and 150 feet (15.2 and 46 meters) above 
mean sea-level within 1 mile (1.61 kilometers) from the spray trial site 
and a maximum elevation of 261 feet (80 meters) occurring at Munden 
Hill, 3 miles (4.8 kilometers) west-southwest from the spray trial site. 
The orchard contains stands of slash pine and Ocala sand pine, some of 
which are stunted and poorly formed because of poor soil and drainage 
conditions. For the purposes of the spray trials, an area on the east 
side of the orchard containing relatively uniform stands of slash and 
Ocala sand pine was selected. Figure 2-1 is a schematic diagram of a 
portion of the selected area. The trees are planted every 15 feet (4.5 
meters) in north-south rows separated by a distance of 30 feet (9.1 
meters). As shown in the figure, slash pine is planted north of the 45- 
foot (14-meter) service road and Ocala sand pine is planted south of the 
road. The locations of the sample trees within which detailed deposition 
measurements were made are indicated by the numbers 1 through 12. Eight 
of the sample trees were located in the slash pine orchard and four were 
located in the more dense Ocala sand pine orchard. A limited survey of 
‘the trees in the vicinity of the sample trees in the slash pine was made 
to locate undersized and missing trees. In the schematic diagram the 
larger open circles represent typical trees and the smaller open circles 
represent undersized or thinly branched trees. The absence of a tree is 
shown by a cross symbol (+). The rows of trees in every case extended 
about 1200 feet (366 meters) to the north and south of the service road 
shown in the diagram. . 


There were 12 aircraft spray trials conducted during the study. As 
shown in Table 2-1, a modified Stearman fixed-wing aircraft was used in 
7 trials and a Hughes 500C helicopter was used in 5 trials. On each 
trial the aircraft flew north and south along the rows or between the 
rows (see Section 3). Except for Trial 14, water containing Nalco-Trol 
and Rhodamine B liquid dye was released on each trial. Nalco-Trol was 
not added to the spray mix on Trial 14. Manganese sulfate was addea to 
the tank mix on Trials 5, 6, 7 and 10 so a mass deposit analysis of 
samples could be performed. 


2.2 Deposition Measurements 


As noted above, deposition measurements were made at the top of the 
canopy, within the canopy and at the ground beneath the canopy. Meas- 
urements at the top of the canopy and within the canopy were made using 
empty aluminum soft-drink cans wrapped on the sides with Kromekote card 
material (4 1/4 inches by 8 3/16 inches) and with a 2 1/2-inch diameter 
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FIGURE 2-1. Schematic diagram of the test area, Withlacoochee State Seed 
Oncnard, Piorida; 19390. 
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TABLE 2-T 


DATES, TIMES AND AIRCRAFT TYPE FOR THE WITHLACOOCHEE STATE 
SEED ORCHARD TRIALS 





Date Time Aircraft 
(1980) | CES BH) Type 


Trial 
Number 





13 Feb i 1554437 | Stearman 








1 
2 15 Feb | 1007: 30 | Stearman 
3 15 Feb | —‘-1506:40 Stearman 
5 16 Feb 1214220 ! Stearman | 
6 16 Feb | Loa? 235 | Stearman 
7 18 Feb | 0942: 24 Stearman 
8 19 Feb | 0734:15 Stearman 
| 10 19 Feb. | :1707:23 | Hughes 500 C 
: ls 20 Feb =| 0748:05 : Hughes 500 C 
12 20 Feb =| ~—-1043:05 | Hughes 00 C 
13 20 Feb : 1756:54 | Hughes 500 C 
14 21 Feb : 0749: 21 Hughes 500 C 
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circle of Kromekote card material taped to the top of the can. At the 
top of the canopy, the cans were suspended from a line stretched east- 
west across the tree rows at the position of the sampling trees. As 
Shown in Figure 2-2 by the dotted lines, the canopy sampling lines 
extended one row east and west of the sample trees and the cans were 
suspended at intervals of 6 feet (1.8 meters) along each line. The 
Kromekote-wrapped cans were also placed within the canopy on the sampling 
trees in the upper crown, mid-crown and in the lower crown at the four 

. cardinal directions (north, east, south and west) for a total of 12 cans 
per sampling tree. In each case, the cans were suspended from the ends 
of branches to simulate pine cones. 


Ground-level measurements were made using 6 11/1l6-by 4 5/16-inch 
Kromekote cards in plastic holders. Ground-sampling lines were placed 
beneath the sampling trees as shown in Figure 2-2 by the solid lines and 
extended two tree rows east and west of the sample trees. Cards were 
placed at intervals of 3 feet (0.9 meters) along these ground-sampling 
lines. Ground measurements using the Kromekote cards in plastic holders 
were also made along the middle of the service road between the slash 
Dine and Qcala sand pine stands at J0-foot ( 3-meter) intervals. The 
service road sampling line began at row 1 in Figure 2-1 and extended 
westward for about 650 feet (198 meters). 


Drift measurements were made using Kromekote cards in plastic 
holders and wrapped soft-drink cans. Depending on the wind direction 
for a given trial, drift samplers were placed along lines to the north, 
east, south or west-of the aircraft spray lines. The east drift line is 
Shown in Figure 2-1. The north and south drift lines were about 1,312 
feet (400 meters) north and south of the service road shown in Figure 2- 
1, just beyond the end of the tree rows in the orchard. The west drift 
line was located about 1,312 feet (400 meters) west of tree row number 
1. If the wind direction for a trial was from, for example, the southwest 
or south, samplers were placed along the north and east drift lines. 
The cards in plastic holders were placed at 100-foot (30.5-meter) intervals 
along the spray drift sampling line and wrapped soft-drink cans were 
placed at 200-foot (6l-meter) intervals on top of 3-foot (l-meter) wood 
Stakes. 


As mentioned above, manganese sulfate was added to the aircraft 
tank mix on Trials 5, 6, 7 and 10 so that atomic absorption techniques 
could be used to obtain direct estimates of mass deposition. Mylar 
sheeting and Kromekote cards were used for subsequent mass analysis at 
selected locations in the ground-sampling network, on all tree soft- 
drink can samplers at 6-foot (1.8 meter) intervals along the canopy 
lines, and-at every other sampling location on the drift sampling lines 
for these trials. 


Not all ground samplers, tree samplers, and top-of-canopy samplers 
were used in all trials because of the time required to put the samplers 
in place. Table 2-2 shows the ground sampling rows and canopy-top 
sampling rows (refer to Figure 2-2 for locations) used in the comparisons 
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FIGURE 2-2. Sampling array for Withlacoochee State Seed Orchard trials, 
Florida, 1930" 


TABLE 2-2 


SAMPLING DATA USED IN THE COMPARISON OF MODEL DEPOSITION 
ESTIMATES AND DEPOSITION MEASUREMENTS 


| Trial Ground Rows eae Rows | 
Number > Tg C-2| C-3 'C-4  C-5 | C-6 | 
2 ; eX. cue XA) paXtahs x 
3 |X Kath | Xo urhl x 
ar et BX “eciXih | wit <ho-X 
Fay ae | 
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TABLE 2-3 


AIRCRAFT DATA 







Speed : Wing or Boom 
: iy Weight No. of 
Aircraft (ms 1) (kg) Se Nozzles* |Nozzles es 





Stearman 40 1406 11.48 D6-46 ag 40 
Hughes 500 C | 11-14 a tey: 8.03 D2-45 (aed 


*Spraying Systems Co. 


of model deposition estimates with deposition measurements presented in 
Section 4. 


2.3 Meteorological Measurements 


Meteorological measurements were made on a 52-foot (16-meter) tower 
located on tree row number 13 at the position shown in Figure 2-1. MRI 
Vectorvanes were mounted at heights of 6.6 and 52.5 feet (2 and 16 
meters) on the meteorological tower. The Vectorvanes measured azimuth 
and vertical wind directions and wind speed which were recorded on 
Esterline Angus strip-chart recorders located at the base of the tower. 
A Climatronics Electronic Weather Station. was used to measure and record 
azimuth wind direction and wind speed at a height of 6.6 feet (2 meters) 
at a position in an open field about 1480 feet (450 meters) west of the 
tree row number 1 and about 66 feet (20 meters) south of the service 
road. Wind direction, wind speed, temperature and relative humidity 
measurements were made to heights of 230 feet (70 meters) during most 
trials using a Tethersonde. The Tethersonde was located about 810 feet 
(247 meters) west of tree row number 1 and 141 feet (43 meters) south of 
the service road in a small open area. 


2-4+-—ircratteData 


Table 2-3 shows the approximate airspeeds and weights of the 
aircraft used in the trials, the aircraft wing or rotor span, and the 
type and number of nozzles used in the spray trials. The spray boom was 
mounted immediately aft of the lower wing on the Stearman. The spray 
boom on the Hughes 500 C helicopter was mounted a few feet below the 
fuselage about midway between the front and rear landing skid struts. 
Table 2-% shows the gallons of material sprayed on each trial, the 
nominal aircraft flight altitude above the top of the canopy, and the 
position of the aircraft swaths flown on each trial with respect to the 
tree row numbers shown in Figure 2-1. 
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BSEGTION Ss 
CBG MODEL INPUTS 


The model inputs required in the CBG model calculations are described 
in this section. The CBG deposition model, used to calculate the deposi- 
tion at the top of the slash and Ocala sand pine canopy, and the CBG 
canopy penetration model, used to calculate deposition within and below 
the canopy, are described in Appendix A. 


Se leecoG Deposition Model Inputs 
Sell Meteorological Inputs 


The meteorological inputs required by the deposition model are 
listed in Table 3-1. The wind speed up, the mean above-canopy wind 
direction 9, and the standard deviations o,{tg} and of used in the 
deposition model calculations were obtained from the MRI Vectorvane 
measurements made at a height of 52.5 feet (16 meters) on the meteor- 
ological tower located at the position shown in Figure 2-1. The Vec- 
torvanes were about 13 feet (4 meters) above the surrounding tops of 
the slash pines. The Esterline Angus strip charts recorded the wind- 
speed and azimuth and vertical wind directions startingl minute before 
the aircraft began spraying the tree rows until 5 minutes after the last 
row was sprayed. On the average, the total data record length for each 
trial was about 12 minutes for the Stearman aircraft trials and about 22 
minutes for the Hughes 500 C trials because of the lower air speed of the 
Hughes 500 C and the need for refilling the spray tanks of this arrcrarec 
during each trial. The wind speed and the azimuth and vertical wind di- 
rections were read from the chart records at 4-second intervals. The 
mean’ wind speed, azimuth wind direction and the standard deviations were 
calculated from the 4-second data records for the first continuous 10- 
minute period of the record after spraying began that did not contain any 
calm wind conditions or missing data. For the eight trials used in the 
comparison of model calculations and measured deposition, the maximum time 
offset for the start of this period was 2.2 minutes after spraying began, 
which occurred on Trial 10. The meteorological deposition model inputs 
for the eight trials derived from the Vectorvane data are shown in Table 
3-2. 


Values of the wind profile power-law exponent p were set equal to 
zero in the deposition model calculations because the aircraft never 
flew at altitudes greater than 19 feet (5.74 meters) above the canopy 
and we believe the data from the top of the meteorological tower best 
represent conditions affecting deposition at the canopy top. The depth 
of the surface mixing layer H. is important in deposition calculations 
when the spray material has relatively small settling velocities and either 
the spray material is released above an inversion, or below a capping 
inversion and deposition is calculated at intermediate and long travel 
distances downwind from the aircraft flight paths. Because the aircraft 
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TABLE 3-1 


METEOROLOGICAL INPUT PARAMETERS 
REQUIRED BY THE DEPOSITION MODEL 


Parameter Parameter Definition 


Standard deviation of the wind azimuth angle (deg) 
for the averaging time es 


Standard deviation of the wind elevation angle (deg) 


Depth of the surface mixing layer’ below a capping 
inversion (m) 


Wind speed at a reference height zp, (m 57!) 
Wind speed profile power-law exponent 


Mean above-canopy wind direction (deg) 





TABLE 3-2 
METEOROLOGICAL INPUTS TQ THE CBG DEPOSITION MODEL 








PPA a) 
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Number i 1) | (deg) | 
al ih ha | TeOZ | 23 : 
9 Coed | ie 20 | 
2 1230 | rhe ! 214 | 
Af 9.9 4.21 | 256 | 
a0) 19.4 6.20 | So 
Ae) 14.6 1276 360 
/, 526 94 305 

6 beae 
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flew just above the canopy and because the gravitational settling veloci- 
ties and aircraft wake settling velocities were significant, the presence 
of an inversion in the vicinity of the canopy would have little effect 

on the deposition. Also, the Tethersonde measurements indicated that 
there were no inversions below a height of 200 feet (60 meters) except 
for Trial 11, when a surface based inversion was present. For these 
reasons a value of 1000 meters was assigned to H_ for all trials to — 
preclude the reflection of drops from capping inversions in the deposi- 
tion calculations. 


Sele Source Model Inputs 


The source model inputs required by the deposition model are defined 
in Table 3-3. The values of Q, H', oo and L used in the deposition 
model calculations are shown in Table 3-4. The source strengths for 
each trial in Table 3-4 were calculated from the expression 


nl Q,\ / 3.78501 + 10° : 

fe a ee (3-1) 

Len gallon ; | 

where | 
Q = Total gallons sprayed (see Table 2-4) | 


as} Number of tree rows sprayed (see Table 2-4) 

In previous studies (Dumbauld, Rafferty and Bjorklund, 1977) the effective 
release height H' has been calculated using an expression accounting for 
the effects of wake turbulence, which cause the spray cloud to sink to a 
height of one-half wing span or rotor diameter (b/2) above the canopy. 
However, because the flight eltitudes were at or below the height (b/2) 

jn the Withlacoochee Spray Trials, the effective release heights H' 

were set equal to the aircraft flight altitude H above the top of the 
canopy. The height (b/2) is also generally used to determine the initial 
source dimensions Dy from the expression 


{xp} = H- (6/2): 
Pema) 


Because the numerator of Equation (3-2) would be negative for the Withla- 
coochee Spray Trials, the following dimensions were arbitrarily assigned 
at the source location (Xp =3 0) \is 
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TABLE 3-3 
SOURCE INPUT PARAMETERS REQUIRED BY THE DEPOSITION MODEL 


’ 





~ Parameter Definition 
Q Source strength (g m4 
al Effective height of the spray cloud release (m) 
Jy {Xp} | Standard deviation of the cloud distribution at the 
| ; 
distance Xp from the source (m) 
L | Length of line source (m) 
‘; Fraction of the total source strength in the ar 
drop-size category 
V5 Gravitational settling velocity for the median drop 
by mass in the j A drop-size category (m s~') 
ns Reflection coefficient for the median drop by mass 
J in the j~ drop-size category 
TABLE 3-4 
INPUT VALUES OF Q, H'; og AND L USED IN 
THE CBG DEPOSITION MODEL CALCULATIONS 
Trial | Ye | be Oy {xa= O} a 
Number | (gbmites)amd | (m) = (m) 
2 ped 95 45 efollilusl 4 4.00 754.6 
3 | 48.28 5.74 4.00 754.6 
5 | 46.26 oc 207 754.6 
6 | BOmS loe 2.6.2 754.6 
7 | 8.36 ia52 2.67 754.6 
10 | 46.15 1.52 1.87 754.6 
1 | 55.74 eee (aay 754.6 
ie | 55.74 We52 1.87 {5556 
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The length of each tree row (2476 feet or 754.6 meters) was assigned to 
the line length L used in Equation (3-1) above. 


The drop-size distributions for use in the deposition model calcu- 
lations were obtained from the drop-size distributions calculated from 
the deposition measurements made on the Kromekote card samplers placed 
on the ground sampling row G-10 located on the service road between the 
slash and Ocala sand pine stands. Briefly, the drop-size distributions 
were obtained by the FS using a Quantimet Image Analyzer and their ASCAS 
(Automatic Spot Counting and Sizing) data program. The image analyzer 
counts the stains produced by the dyed spray drops and measures drop- 
stain diameters on a representative area of the sampling card in up to 
16 pre-selected size categories. The ASCAS program uses the counts in 
the 16 size categories, a mathematical equation relating stain diameters 
to actual drop diameter, and the drop density (1g cm73, in this case) 
to produce the mass distribution measured on a sampling card or group of 
cards. Table 3-5 contains the values of F5> V; and Yj assigned to the 
drop-size categories measured for Trial 2. Similar information for al] 
eight trials used in this-study is contained in Appendix B. The values 
of f; were obtained directly from drop-size distribution information 
produced by the ASCAS program. Values of the gravitional settling 
velocity Vs; were calculated for drops of unit density and the mean drop 
diameters tor each category using the technique described by McDonald 
(1960). The values of the reflection coefficient Y; shown in the table 
are based on the relationship between gravitational settling v.locities 
and reflection coefficients postulated by Dumbauld, Rafferty and Cramer 
(1976). A value for y; of 0 means the drops in this category are deposited 
at the ground or, for deposition calculations at the canopy top, that 
all drops in the category enter the top of the canopy. A value of botog 
y.means that all drops are reflected at the canopy top or at the ground. 

e value for Y; of 0.59 in Table 3-5 for the first drop-size category 
means that 59 percent of the drops reaching the canopy top are reflected 
and 41 percent enter the canopy. 


3.2 CBG Canopy Penetration Model Inputs 
th pal Meteorological Model Inputs 


The mean wind speeds in quarters of the forest canopy are the only 
meteorological inputs required by the canopy penetration model. Although 
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, TABLE 3-5 
DROP-SIZE DISTRIBUTION, SETTLING VELOCITIES V. AND 


REFLECTION COEFFICIENTS Yj FOR TRIAL 2 OF THE 
WITHLACOOCHEE SPRAY TRIALS 


Drop-Size Mean Drop | 3 V 





Category ; Diameter 5} J ‘5 
j + Gum) J (ns!) 
1 | 38 .0 0.001 0.0629 0.59 
2 i 7st 0.009 0.216 0.21 
3 | 118.0 0.02 0.402 0 
4 | 153.0 0.03 0.595 0 
5 i 185.0 0.04 0.741 a: 
6 e522980 0.10 0.926 0 
7 P85 78¥0 0.10 1.20 0 
g | 327.0 0.10 1.44 0 
9 406.0 0.20 1.76 0 
10 | 486.0 0.10 2.07 0 
11 | 546.0 0.10 2.36 0 
12 623.0 0.10 2173 0 
13 705.0 0.04 3.11 0 
14 792.0 0.03 3.44 0 
15 933.0 0.02 3.96 0 
16 jal Aah 0.01 4.93 0 
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as mentioned in Section 2, a Yectorvane was located at a height of 6.6 
feet (2 meters) on the meteorological tower, the wind speeds within the 
canopy were usually so light and variable that the accuracy of the mea- 
Surements was questionable. For this reason, we used the normalized 
wind profile given by Fritschen, et al. (1970) for an isolated conifer 
stand of similar stem density with no understory to define the below- 
canopy wind profile. The below-canopy wind speeds in Table 3-6 were 
obtained by assuming that the wind speed measured at the top of the | 
meteorological tower was equivalent to the value indicated by the Fritschen 
profile at the canopy top. The wind speeds in Table 3-6 were used in 
calculating the deposition below the canopy described in Section 4. 


Shae Canopy Inputs 


The canopy penetration model requires the canopy data inputs and 
related information described in Table 3-7. The values used in the 
Canopy penetration model calculations for this study for Hc, Dt, 

PRPEN and M are given in Table 3-8. Tree widths are given in Table 3-9. 
The average tree widths and heights were determined from measurements of 
the sample trees made during the trials. The values of PRPEN were 

based on estimates provided by Grim and Barry (1975) for conifer trees 
with moderate (their foliage type III for the slash pine) and heavy 
(their foliage type IV for the Ocala sand pine) foliage densities. The 
impaction collection efficiencies for the canopy penetration calculations 
were calculated from the following empirical relationship, recommended 
by Grim and Barry (1975) and attributed to Sell: 


Ae} eee Ge | 
fo Se Ne Fr) ee 2 ere a} 
us : 
E. = ; (3-4) | 
J | | 
1 : E Doel 
where 
u = impaction velocity in meters per second 
q, = drop diameter for the jth category in micrometers 
S = diameter in centimeters of the element on which the drop impacts 


Grim and Barry obtained a good fit to their deposition data using a value 
for s of 13 centimeters. We have also used this value of s in model 
calculations made for the Rennic Creek Trials (see Dumbauld, Rafferty 

and Bjorklund, 1977). The impaction velocity u in Equation (3-4) was 

set equal to the mean wind speed in the forest canopy obtained by averaging 
the four values for each trial shown in Table 3-6. Values of E; for 

each drop-size category are shown for each trial in the tables contained 

in Appendix B. ; 
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TABLE 3-6 


MEAN. WIND SPEEDS IN METERS PER SECOND USED IN THE CANOPY 
PENETRATION MODEL CALCULATIONS 


Canopy Height Interval (k) 





Trial 
Number ] 2 3 
(Base) 
| 2 0.95 0.84 0.82 Q3 
| 3 ¥-02 0.90 0.88 We 
| 5 2.93 2050 ZeeV “at 
| 6 3292 3.45 checy! 3: 
7 B74 5.06 ~ 4,94 Be 
10 1.64 1.44 1.41 es 
1] 0.87 Oo77 Oe75 0: 
i ALS, 1256 1.52 1 
TABLE 3-7 
CANOPY RELATED INPUTS REQUIRED BY THE 
CANOPY PENETRATION MODEL 
Parameter Description | 
H Tree height (m) 
0. Tree density (stems per acre) 
PRPEN Probability of penetration 
M Total number of drops to be passed along the 
trajectory 
W. Tree width at one-meter height intervals (m) 
Ee Impaction collection efficiency for the fia 


drop-size category 





ie) 


TABLE 3-8 


VALUES OF CANOPY RELATED INPUTS USED IN THE 
CANOPY PENETRATION CALCULATIONS 





Parameter 









Slash Pine Ocala Sand Pine 


| 

: D, (stems/acre) 96.8 | 
| PRPEN 0.13 
| M (drops ) 1000 : 
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TABLE 3-9 


TREE WIDTHS FOR THE SLASH 
AND OCALA SAND PINES 


Tree Width (m) 


Slash Pine Ocala Sand Pine 


1 
2 
3 
4 
5 
6 
7 
8 
9 
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: SECTION 4 
RESULTS OF THE CALCULATIONS 


The model inputs described in Section 3 were used in the line-- 
source deposition model presented in Appendix A to calculate deposition 
at the canopy top along sampling rows C-1 through C-6 (see Figure 2-2). 
Deposition calculations along the ground-level sampling rows G-1 through 
G-9 were made using the canopy-penetration model described in Appendix 
A. The results of these calculations and comparisons of these results 
with observed deposition are described below. 


4.1 Graphs of Observed and Model-Calculated Deposition 


Graphs of observed and model-calculated deposition at the top of 
the slash pine canopy and on the ground beneath the canopy for all 
trials are presentéd in Appendix C. Figure 4-1 shows the observed and 
model deposition at the canopy top for Trial 2 in which the Stearman 
aircraft flew at a height of 19 feet (5.74 meters) above the canopy. In 
Figure 4-1, downline distance is measured from west to east along the 
sampling rows C-1 through C-4. The observed deposition plotted in the figure 
is the arithmetic average deposition at the sampler positions located at 
the same downline distances. For example, the observed deposition — 
values at samplers C-1-1, C-2-1, C-3-1 and C-4-1 were summed and divided 
by 4 to obtain the deposition plotted at the downline distance of Ose 
The model deposition values shown in Figure 4-1 are-also average deposi- 
tion values at each downline distance. 


Inspection of Figure 4-1 shows that the peak depositions along the 
sampling rows predicted by the model are greater than the observed peaks 
and the model depositions in the troughs are much smaller than the ob- 
served depositions in the troughs. The peaks in the model and observed 
depositon pattern reflect the maximum deposition at the top of the can- 
opy from individual swaths of the Stearman aircraft. Our experience in 
modeling aircraft spray deposition patterns indicates-that the ratio 
between the model depositions for the peaks and troughs is principally 
controlled by the value assigned to the source dimension og at the 
height of cloud stabilization above the canopy. In previous comparisons 
of observed and model-calculated deposition patterns for aircraft flying 
at several wingspans above the top of the canopy (Dumbauld, Rafferty and 
‘Cramer, 1976; Dumbauld, Rafferty and Bjorklund, 1977), the use of Equation 
(3-2) to determine og yielded satisfactory agreement between the 
observed and model-calculated peaks and troughs in the deposition pat- 
tern. As noted in Section 3.1.2 above, Equation (3-2) is based on aircraft 
wake effects and could not be used for the Withlacoochee Spray Trials be- 
cause the aircraft flight altitudes were less than half a wing span or 
rotor diameter above the canopy. Therefore we used an arbitrary expression 
for ag given by Equation (3-3). The effects of doubling the value of 
Tg assigned on the basis of Equation (3-3) for the Trial 2 model calcula- 
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tions are illustrated in Figure 4-2. Comparison of Figures 4-1 and 4-2 
shows that doubling the value of og produces a model deposition pattern 
more nearly reflecting the peak-to- trough ratio in the observed deposi- 
tion pattern. The doubling of og also increases the average calculated 
depos Loan level along the downline distance from 3.61 x 103 to°4:75 x 

103 mg m-2, as compared to the average observed deposition level of 5.12 
x 103 mg m-2. Further comparisons of average observed deposition with 
model average deposition along the downline distances are-provided in 
Section 4.2. 


Figure 4-3 shows calculated and observed deposition at ground level 
beneath the slash pine canopy for Trial 2. The averaging techniques are 
the same as those used to obtain the results shown in Figure 4-1 for 
observed deposition and were applied to rows G-l, G-3, G-4 and G-6. 
Figure 4-3 shows that the peak-to-trough ratios of the model calculations, 
also made using og based on Equation (3-3), are smaller than the 
corresponding ratios for the top of the canopy shown in Figure 4-1 but 
larger than the ratios indicated by the observed deposition. Also, the 
average model deposition along the downline distance is in better 
agreement with the observed average ground-level deposition than for the 
top of the canopy. 


After Trial 2, the Stearman flew at an altitude of only 5 feet (1.5 
meters) above the canopy. The Hughes helicopter flew at an altitude of — 
5 feet (1.4 meters) above the canopy on all trials. Consequently, 
Equation (3-3) yields smaller values of ao, than for Trial 2. Figure 
4-4 shows observed and calculated deposition- patterns at the top of the 
canopy for Trial 3. Trial 3 was conducted with the Stearman under 
meteorological conditions similar to Trial 2, except that vertical 
turbulence (indicated by o¢ in.Table 3-2) was greater. The use of 
larger values of o¢ in the model calculations would, if atmospheric 
dispersion dominated the deposition process at the top of the canopy, 
act to decrease the peak-to-trough ratios in the model-calculated deposi- 
tion patcern. However, Figure 4-4 shows that the calculated deposition 
pattern still exhibits much greater peak-to-trough ratios than the 
observed pattern, indicating that the value assigned to og is the 
dominant factor. The figures in Appendix C for deposition at the canopy 
top for the helicopter trials (Trials 10, 11 and 13) also show that we 
have underestimated the value of og in the model calculations. The 
peak-to-trough ratios for the model patterns for Trials 6 and 7 in 
Appendix C agree better with the observed ratios, possibly because the 
mean wind speeds for these trials were greater. 


The. value of the initial source dimension og is principally 
controlled by the aircraft wake effects. The wake effects_of the rela- 
tively slow flying helicopter (air speeds of 11 to 14 m s~*) are much 
greater than those for the Stearman. We believe the discrepancy between 
the peak-to-trough ratios of the observed and calculated patterns is due 
to our inability to define 9% as a function of the fixed-wing and 
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1980, using an initial source dimension O, twice as 
large as that used in the model calculations snown in 
Figure 4-1. 
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helicopter wake effects for these low-level flights. The value of do 
used in the model calculations also affects. the average deposition 
levels along the downline distances, but to a lesser degree. 


4.2 Ratios of Model to Observed Average Deposition Levels 


A measure of the overall model performance in calculating deposition 
at the top of the slash and Ocala pine canopies, and at ground level 
beneath the trees, is obtained by forming the ratio of model to observed 
average deposition along the sampling lines. Both the model and observed 
average deposition levels were obtained by summing deposition levels 
over all samplers and dividing by the total number of samplers. The 
results are presented in Table 4-1. The value of 0.80. shown for the 
ratio at the slash pine canopy top for Trial 2 is the ratio of the mode] 
average deposition level along the downline distance (5.776.x 103 mg m-2) 
to the observed deposition level (7.219 x 103 mg m-2) shown in Figure 4-1. 


Table 4-1 shows the observed average deposition levels are greater 
than model depositon levels at the slash pine canopy top on all trials 
and at the Ocala sand pine canopy top on all trials except Trials 1 and 2. 
Thus, there appears to bea bias in the model towards underprediction of the 
deposition at the canopy top. We believe that this is caused by the failure 
of the model to account for aircraft wake effects. Also, the ratios for 
trials conducted under similar conditions fluctuate over a large range. 
For example, the ratio in Table 4-1 for Trial. 5 at the top of slash pine 
canopy is 0.25 and for Trial 6, conducted under nearly identical conditions, 
the ratio is 0.75. We have compared the average observed deposition 
with the average deposition that should have occurred if all the material 
released by the aircraft deposited uniformly within the spray grid area, 
assuming no spray drift took place. The results are shown in Table 4-2, 
where the average mean deposition for total recovery TR is calculated 
From the expression 


1000 0 | 
asia (4-1) 


—_— 


where 


Q = source strength in gm (see Table 3-4) 


SW = swath width or distance between tree rows for these 
teialss( 9. im) 


The ratios of observed deposition levels to the calculated levels based 
on total recovery clearly show that, for most trials, the observed depo- 
Sition levels obtained from the Kromekote card samplers indicate more 
material was deposited than appears reasonable. Difficulties in the 
analysis of the Kromekote card data were reported by FMP/MAG due to 
neavy drop density on the cards and, in some cases, spreading of drop 
stains because the cards were dampened by dew or mist. Under these 
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TABLE 4-1 


RATIOS OF MODEL TO OBSERVED AVERAGE DEPOSITION ALONG THE 
DOWNWIND SAMPLING LINES* 


Trial Slash Pine Canopy Ocala Sand Pine 


Top Ground Ground 


Number 


Geometric 
Mean Ratio 





* Based on mass deposition 
**x No measured data 
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TABLE 4-2 


AVERAGE OBSERVED DEPOSITION AND AVERAGE DEPOSITION BASED ON 
TOTAL RECOVERY AT THE CANOPY TOP 


Average Observed Deposition Total Ratio 
SIME TR 
Slash Pine |Ocala Sand Pine (mg m Slash Pine jOcala Sand Pine 


(mg m-¢) Observed/Total Recovery 


T0229 exe 10 
Rye 9% 
.490 x 
* 

peed 

. 386 

Mes 2 

. 296 





* 
No samplers were placed in the Ocala sand pine canopy on this trial. 


30 





circumstances, the Quantimet Image Analyzer may have made innacurate 
measurements in counting and sizing drop stains, due to the high drop 
densities, that lead to both overestimation of mass deposition and 
variability in the ratios. 


Since it is unrealistic to expect the model to predict deposition 
levels at the canopy top greater than those based on total recovery, we 
nave constructed adjusted ratios of model average deposition levels to 
"observed" deposition levels. The adjusted ratios were obtained using 
the deposition levels based on total recovery given in Table 4-2 as 
observed deposition levels when the ratios in Table 4-2 indicated more 
material was deposited than appears reasonable (ratios greater than 
unity in Table 4-2). The results of adjusting the ratios are shown in 
Table 4-3. We believe the ratios given in Table 4-3 more accurately 
reflect model performance in predicting average deposition levels. 
Furthermore, we would expect the observed average deposition levels at 
the top of the canopy to be less than those based on total recovery. If 
the observed deposition levels were less than those based on total 
recovery, the geometric mean ratios in Table 4-3 would be increased. 


Drop densities on the cards beneath trees on the ground sampling 
rows were generally less than those at the top of the canopy, although 
between rows there were some very high drop densities and moisture on 
the cards was a problem in some trials. If the counting and sizing of 
drops on these cards were more accurate and thé canopy penetration model 
performs well, we would expect the ratios for the ground sampling rows 
in Table 4-1 to compare favorably with the adjusted ratios in’Table 4-3 
for the top of the canopy. Comparison of these ratios in the two tables 
indicates that the ratios for the slash pine (0.76 and 0.75, respectively) 
are comparable. However, the geometric mean value for the ground sampling 
rows beneath the Ocala sand pine (0.32) is much less than the geometric 
mean ratio for the Ocala sand pine canopy (0.80) in Table 4-3. This 
may indicate, for example, that the value assigned to the probability 
of penetration (PRPEN = 0.13) for the Ocala sand pine canopy was too 
small while the value assigned to the slash pine canopy (PRPEN = 0.38) 
was more realistic. However, the large trial-to-trial variation of the 
ratios in Table 4-1 for the ground sampling rows beneath the Ocala sand 
pine canopy most likely indicate that other factors are involved. 


As noted in previous sections, Mylar samplers were used on Trials 
5, 7 and 10 and an atomic absorption technique used to analyze for 
manganese salts in the washed samples. Deposition estimates based on 
these measurements were available for six locations along each ground 
and canopy-top sampling line for these trials. The canopy-top samplers 
(soft drink cans) had Mylar on the top of the can as well as around the 
sides of the can. The Mylar from the top and sides of the can was 
washed and analyzed simultaneously. Since the deposition model calcula- 
tions are only strictly comparable with deposition measurements from the 
top of the can, comparison with the observed deposition measurements 
based on this analysis is not possible. The ground Mylar samplers were 
placed in Kromekote card holders and the results of the atomic absorbtion 
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TABLE 4-3 


ADJUSTED RATIOS OF MODEL TO OBSERVED AVERAGE DEPOSITION 
LEVELS AT THE CANOPY TOP 


Ocala Sand Pine 
Canopy 


Slash Pine: 
Canopy 









Geometric 
Mean Ratio 





No samplers were placed in the Ocala sand pine canopy on this trial. 


oye 


analysis can be directly compared with the results from the ASCAS 
analysis of the Kromekote cards and with the mode] calculations. Table 
4-4 shows average observed deposition values obtained by arithmetically 
averaging the deposition values available from all Kromekote cards and 
Mylar samplers beneath the slash and Ocala sand pine canopies. The 
Mylar samplers were placed next to the Kromekote card samplers at 
positions 30 through 35 along the sampling lines (approximately in the 
center of the sampling line in the vicinity of the sampling trees) and 
Spanned a distance of 4.6 meters along the sampling line. The average 
observed values for the Mylar and Kromekote samplers beneath the slash 
pine canopy agree reasonably well for Trials 5 and 7 and beneath the 
Ocala sand pine canopy for only Trial 7. There is very poor agreement 
between the two measurement techniques on Trial 10. In fact, the Mylar 
Samplers show paatae?) deposition levels greater than expected from total 
recovery (5047 mg m Since the Mylar samples spanned a relatively 
Short distance along ic sampling line, spatial positioning of the 
model-calculated peaks and troughs with respect to the line and the peak- 
to-trough ratios of the model pattern can greatly affect the model 
comparison with the measurements. For example, in Trial 10 for the 

Slash pine canopy, the model pattern has sharp peaks and troughs because 
of a failure adequately to account for the wake effects of the helicopter. 
As a result, the model shows the Mylar samplers to be closer to a trough 
in the calculated pattern while the measurements indicate that they are 
closer to a peak in the pattern. Therefore, the model greatly under- 
predicts the observed calculations indicated by either of the measure- 
ment techniques. For the samplers beneath the slash pine canopy on 
Trials 5 and 7 where the two measurement techniques agree and the model 
pattern does not show extreme ranges in deposition between the peaks and 
troughs, the model calculations and the measurements are in closer 
agreement. 


The results of the calculations described above are summarized in 
Section 5. 
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TABLE 4-4 


AVERAGE OBSERVED (MYLAR AND KROMEKOTE SAMPLERS) 
AND MODEL DEPOSITION LEVELS (mg m-¢) FOR 
GROUND SAMPLERS BENEATH THE CANOPY* 


Trial Gano Average Observed Average Observed Model 
Number PY (Mylar Sampler) (Kromekote Sampler) Average 
986 | 928 


278 









* The values given in the Table for Kromekote samplers and the model 
calculations are based on deposition values for the positions where 
both Mylar and Kromekote samplers were located. 
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SECTION 5 
SUMMARY AND RECOMMENDATIONS 


5.1 Summary of Results 


The CBG model computer program has been used to calculate deposition 
patterns for comparison with measurements made at the canopy top and on 
the ground beneath the canopy in the Withlacoochee spray trials. The 
deposition patterns at the canopy top discussed in Section 4.1 show that 
the peaks and troughs of model deposition patterns for the fixed-wing 
aircraft are exaggerated compared to the measured patterns, especially 
under low wind-speed conditions. The model peaks and troughs for the 
helicopter deposition patterns are also greatly exaggerated compared to 
the observed patterns. At ground level, the model deposition patterns 
for the fixed-wing aircraft more nearly approximate the observed deposi- 
tion patterns, while the model helicopter patterns do not. We believe 
the failure of the model calculations to match the shape of the observed 
patterns is caused by aircraft wake effects that are not adequately 
accounted for in the model. In almost all the Withlacoocnee trials, 
both aircraft flew at a height of about 1.5 meters above the canopy. 

The wake effects expressions in the CBG model assume that the minimum 
aircraft altitude above the top of the canopy is greater than one-half 
the wing span (11.48/2meters) or rotor diameter (8.03/2 meters). In all 
previous tests of the performance of the CBG model in calculating 
deposition patterns from fixed-wing spray aircraft and helicopters 
(Dumbauld, Rafferty and Cramer, 1976; Dumbauld, Rafferty and Bjorklund, 
1977) the above-canopy flight altitudes were larger than the minimum ~ 
altitude assumed in the wake effects expressions and the shapes of the 
model and observed deposition patterns were similar. 


Table 5-1 shows geometric mean ratios of model to observed average 
mass deposition at the canopy top of the slash and Ocala sand pines and 
at the ground beneath the canopies. ‘According to these ratios, the CBG 
model program underestimates the deposition at the top of the slash pine 
canopy by more than a_ factor of 2, and underestimates the deposition at 
the ground by a factor of 1.3. The model. also underestimates the average 
deposition at the canopy top of the Ocala sand pine by a factor of 2 and 
at the ground by a factor of 3. However, the observed average depositions 
reflected in these ratios are larger than expected if all the material 
emitted by the spray aircraft entered the canopy. We believe this anomalous 
result is explained in part by limitations of the Quantimet drop counting 
and sizing analysis. High density of drops on the sampling cards and 
the presence of moisture on the cards during some trials resulted in an 
overestimation of the mass deposition on the cards. This was due to in- 
ability of the Quantimet to distinguish single drops from drops which 
touch. Also, dampness causes excessive spreading on Kromekote cards and 
thus an overestimate of drop size and spray volume. For these reasons, 
we adjusted the average deposition levels at the canopy top so that they 
did not exceed the values expected from total recovery based on the ap- 
plication rate Q (see Table 3-4). Table 5-2 shows adjusted geometric 


TABLE 5-1 


GEOMETRIC MEAN RATIOS OF AVERAGE MODEL TO OBSERVED 
MASS DEPOSITION FOR THE WITHLACOOCHEE SPRAY TRIALS 


Canopy Top Ground 


Fixed Wing He1icopter| Combined |Fixed Wing|Helicopter |Combined 


Slash 
Pine 


Ocala 
Sand Pine 








TABLE 5-2 


ADJUSTED* GEOMETRIC MEAN RATIOS OF MODEL TO 
OBSERVED MASS DEPOSITION AT THE CANOPY TOP 


ee Dea Fixed Wing Helicopter Combined 
Slash 
Bh ne 0.76 Biase 
Ocala ; 
Sand Pine Q.95 0.80 
Combined 0.84 Dial y 





* 
Adjusted to account for possible errors in counting and sizing drops 
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mean ratios of model to observed average mass deposition using the ad- 
justed levels. The model calculations underestimate the observed average 


deposition at the canopy top by about a factor of 1.3 after this adjust- 
ment was made. 


Since the adjusted geometric mean ratio for the slash pine canopy 
top in Table 5-2 (0.75) agrees with the ratio for ground level beneath the 
canopy in Table 5-1, it appears that model inputs to tne canopy penetra- 
tion model (Grim and Barry, 1975) in the CBG computer program have been 
chosen correctly. On the other hand, the difference of the ratios for 
the Ocala sand pine at the top of the canopy and at the ground indicates 
that the inputs to the model for Ocala sand pine require further adjustment. 


Manganese salts were added to the spray material during 3 trials as 
a chemical tracer to permit a subsequent direct analysis of mass deposi- 
tion using atomic absorption techniques. Comparison of observed average 
deposition at the canopy top and beneath the canopy based on this analysis 
technique with the average deposition obtained from the drop counting 
and sizing technique show poor agreement between the two analysis methods. 
Only a limited number of the Mylar samplers were used to collect spray 
deposition for the atomic absorption analysis and this technique was 
used only in 3 trials. For this reason there were insufficient deposition 
data based on this analysis technique for adequate assessment of the 
model performance. 


In spite of the expected difficulties in applying the CBG model for 
use with low-flying aircraft, we believe the results of this study demon- 
strate the feasibility of adapting the CBG model for applications in 
planning and conducting aerial application programs in seed orchards. 


5.2 Recommendations 


As explained above, serious problems were encountered in adequately 
accounting for the wake effects of low-flying aircraft ard in the counting 
and sizing of drops under the experimental conscvraints of the Withlacoo- 
chee Spray Trials. The results clearly show that improvement in the 
wake model, used to obtain source dimensions, is required to adequately 
describe the shapes of deposition patterns when the aircraft flight 
altitude is in close proximity to the top of the canopy. Until this 
improvement is made, the utility of the model in defining optimum spray 
aircraft altitudes and swath widths for orchard spray operations is 
limited. The estimation of spray drift, specification of buffer zones, 
and effective deposition levels within the canopy depends in part on the 
ability to calculate accurately the deposition at the top of the canopy. 
To properly assess the utility of the CBG model computer program in 
calculating average deposition requires more quantitative deposition 
measurements than those made during the Withlacoochee trials. 


Because of funding and time limitations, we have been unable to 


utilize all the data available from the Withlacoochee trials to improve 
model performance. For example, larger values for the model input 


a 


PRPEN describing the probability of a drop penetrating a tree may result 
in more accurate estimates of ground-level deposition and canopy retention 
of drops for trees similar to the Ocala sand pine, which had denser fo- 
liage than the slash pine. 


We make the following recommendations for future work: 


8 -fforts should be made to improve the wake model portion 
of the CBG program concept for low-flying aircraft 


@ The manganese salt tracer technique should be further - 
assessed as a means of measuring mass deposition in 
cases where Kromekote cards are expected to show high 
drop densities 


a Further analysis of the Withlacoochee trial data should 
be performed to improve model parameteric inputs and 
reassess model performance where practical 
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APPENDIX A 
DESCRIPTION OF THE CBG DEPOSITION AND 
CANOPY PENETRATION MODELS 


A.1 CBG Deposition Model 


In the model, the axis of the spray cloud is assumed to intersect 
the ground at a downwind distance from the source that is proportional 
to the_product of the release height H and the mean cloud transport 
speed u divided by the drop settling velocity Vj for the jth drop-size 
category. The inclination of the cloud axis from a ees for a 
given drop-size category j is thus equal to tan™ It is assumed 
that drops dispersed upwards by turbulence are meee esa down- 
wards at the top of the surface mixing layer, but the fraction of drops 
reflected at the ground is a variable input parameter for each j 
category. The model uses the Cartesian coordinate system shown in 
Figure A-l for a line source of length L at a height H and a calculation 
point at R(¢€, 5, z). Deposition, expressed in units of mass per unit 
area, at the point (¢, 95, 0) is given by the expression 
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(Equation (A-1) continued on following page.) 
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FIGURE A-1. Schematic diagram showing the line source geometry with 
respect to a calculation point at R (ec, §, z) and wind 
direction 4. 
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(Equation (A-1) Continuéd. ) 
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The following parameters used in the preceding equations are based 
On meteorological measurements or inferred from meteorological measurements: 


Tn = standard deviation of the wind azimuth angle in radians 
kK = constant relating om and or 
rs (A-19° 
Of = standard deviation of the wind elevation angle in radians 
Hn = depth of the surface mixing layer below a capping inversion 
U = mean transport wind speed 
U ea - zl 
Rak 2 1 saietebar 
i p R 
aaherite)aicee) 
; (A-20° 
Up sel EO IUS 
Up = mean wind at the reference height z, 
p = wind power-law coefficient 


A=5 


Zo = effective upper bound of the cloud 
aber eo kxe tx ft ra eee a 
; ory Ve 2 m 
= (A-21) 
Hn ; 25 S a 
Z4 = effective lower bound of the cloud 
KXe 
Hepeanee eel Ga F *y) > 2. Z 
= (A=22 
2 5 2 < 2 


angle between a line perpendicular to the line source 
and the mean wind direction (see Figure A-1) 


The following parameters are source inputs requ red for use in the 
model: 


Q = total source strength emitted along the length L of the 
line source 


H = release height 


V. = gravitational settling velocity for the median drop by 
mass in the jth drop-size category 


f. = fraction of the total source strength in the yun drop-size 
category 


A-6 








y. = reflection coefficient for the median drop by mass in the 


J jth drop-size category 
Oh standard deviation of the cloud distribution at the dis- 
tance x 
R 
Xp = distance from the line source to cloud stabilization 
L = length of the line souce 


Despite the apparent complexity of the model, solution times are 
relatively short because the summation for multiple reflection from i=0 
to infinity usually can be terminated after less than 10 passes (by 
checking to see if succeeding passes result in significant increases in 
the calculated deposition). Also, as shown by the model equations, the 
solution fails when the wind is exactly perpendicular to the line source 
(6=0 degrees) or when the wind direction is exactly parallel to the line 
source (8=90 degrees). Exact solutions for these angles are easily 
derived. However, it is simpler to have the computer program change 
these angles by a small fraction (0.01 degrees) and then calculate the 
deposition because there is no significant loss in accuracy. 


A.2 CBG Canopy Penetration Model 


The model is based on a Monte Carlo technique where a large number 
of drops in each size category are passed along a trajectory through 
a simulated forest with trees assigned to equal areas according to the 
density (stems per acre) in the forest being simulated. The drop trajec- 
tory is a function of the gravitational settling velocity Vj and the 
mean wid speed at various levels within the forest canopy. As a drop 
proceeds along the trajectory, each tree is randomly displaced within 
the assigned area in the plane of the horizon and calculations are made 
to determine if the drop intersects the tree envelope and, if an inter- 
section occurs, whether the drop strikes a tree element. When the drop 
strikes a tree element, a tally is recorded for the height interval 
within the canopy where the “hit" occurs and for all greater height 
intervals. Drops proceed along the trajectory until a hit occurs or 
until the trajectory intersects the ground. After the specified total 
number of drops in the size category have passed along the trajectory, 
the tally number within each height interval is divided by the total 
number of drops to obtain the percentage of drops reaching the given 
height interval. The total number of drops passed along the trajectory 
required to achieve a stabilized solution (percentage penetration) is a 
function of the steepness of the trajectory, with more drops being 
required for size categories with large settling velocities. 
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Figure A-2 is a schematic diagram showing an example drop tra- 


jectory and forest construct. 


intervals is defined by the following piecewise linear function: 
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U 
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0 3 k = 0 
0.25(k-1)H 
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angle, measured in radians, defining the 
declination of the trajectory with respect 
to the plane of the horizon 





-] Vie 
tan = J : 'o < 1.4 radians 
: Uosk > 
1.4 radians ; o> 1.4 radians 


height of the forest canopy 


mean wind speed in the yeh 


within the canopy 


height interval 


The number of trees placed along the trajectory path in the simulated 
forest is given by the expression 


N = Xq/OX, 
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The drop trajectory in each of k height 


(A- 23) 


(A-24) 
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(A-26) 
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where 


Xa = maximum horizontal travel distance in meters 
of the drop within the forest canopy 


AX, = alongwind tree spacing within the simulated 
forest 
= 63.8//D, 
D, = tree density in units of stems per acre 


As shown in Figure A-2, each tree stem is given the following 
location along the trajectory: 








X, x (N, - n) Ax, + (R_-0.5) Ax, 5  n=l,2,...-N, 
Ye = (R' - 0.5) AX, 
where 
n = tree number 


and R and R' are uniform random numbers between 0 and 1. The possi- 
bility of the drop trajectory intersecting the tree envelope is deter- 
mined by comparing the distance yt with the radius of the tree envelope 
(W+/2) at the height Zp where the trajectory passes through the distance 
Xt. The tree widths W, are specified by the program user at one-meter 
intervals and the program calculates the radius (W+/2) at Zp by linear 
interpolation., If yt is greater than (W+/2), no intersection occurs and 
the computer program proceeds to the next tree along the trajectory. 


If the drop trajectory intersects the tree envelope, the program 


calculates the probability that the drop impacts on a tree element from 
the expression 


Bhs Bal PRPENC ) 
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where 


Ee = pe pecen Oe efficiency of the tree element for the 
J jth drop-size category 
PRPEN = probability of penetration for the population of 
drops and for a horizontal trajectory through the 
tree 
(pa path length correction factor for a non-horizontal 
trajectory 
pal a : ek 
COS , Wa 
Cr Ss (A-32) 
ue a ls 
Wo Ws 
i = maximum width of the tree envelope 


A particular drop of the population of drops'is assumed to intersect the 
tree element when the value of Pj from Equation (A-31) is greater than a 
uniform random number R" between 0 and 1. Each tree is divided into ten 
height-class intervals, and an intersection with a tree element is. 
recorded as a "hit" in the class interval in which Zp Occurs and in 
every higher class interval. 


The process described above is repeated for every drop passed along 
the trajectory and the final percentage of material penetrating- to a 
given height interval determined by dividing the number of recorded hits 
in the height interval by the total number of drops in all jth size 
categories passed along the trajectory. 


The inputs required by the canopy-penetration model are: 


PRPEN = probability of penetration 


Uae mean wind speed in the Kin height interval within the 
i. forest canopy 
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; = gravitational settling velocity in meters Be yr 
second for the median drop by mass in the jth 
drop-size category 


th 


E, = impaction collection efficiency for the j~ drop- 
size category 
D. = tree density in stems per acre 
We = tree height in meters 
W. Sd tree width at one-meter height intervals 
M = total number of drops to be passed along the 
trajectory 


The computer program also permits the user to simulate a multi-storied 
canopy of up to three tree heights, with a different value of PRPEN 
possible for each tree or story height. 


The output from the canopy penetration model includes, in addition 
to the percentage of. material in the jth drop size category penetrating 
to each of ten levels within the canopy, the maximum horizontal travel 
distance of the drop within the forest canopy (xq). When the deposition 
model given by Equation (A-1) is used in conjunction with the canopy 
penetration model to calculate ground-level deposition within the 
canopy, the distance from the line source to the target is adjusted so 
that the deposition model calculates the deposition at the top of the 
forest canopy at a distance xq upwind from the target. The deposition 
at the target is then determined by multiplying the deposition predicted 
from Equation (A-1) by the percentage of material reaching the ground 
from the canopy penetration model. 
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APPENDIX B 


This appendix contains tables of drop-size distributions and associated 
gravitational settling velocities, reflection coefficients and impaction 
efficiencies for Trials 2, 3, 5, 6, 7, 10, 11 and 13 of the Withlacoochee 
Spray Trials. 
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TABLE B-1 


DROP-SIZE DISTRIBUTION, SETTLING VELOCITIES Yj, REFLECTION 
COEFFICIENTS Y; AND IMPACTION EFFICIENCIES Ej FOR 
TRIAL 2 OF THE WITHLACOOCHEE SPRAY -TRIALS, 
FLORIDA, 1980 





Ei Drop-Size Mean Drop ; 
i Category Diameter Y. E.. 
| j ( um) J : 
1 38.0 0.001 0.063 0.59 0.027 
2 78.3 0.009 0.216 0.21 0.117 
3 118 0.02 0.402 0 0.265 
4 153 0.03 0.595 0 | 0.445 
5 185 0.04 0.741 0 | 0.650 
6 228 0.10 0.926 ao 0.988 
7 278 0.10 1.20 0 1 
8 327 0.10 1.44 0 1 
9 406 0.20 1.76 0 1 
10 486 0. 105° ZnO] 0 1 
11 546 0.10 2.36 0 1 
12 623 0.10 oa3 0 1 
13 705 0.04 cee 0 1 
14 792 0.03 3.44 0 1 
15 - 933 0.02 3.96 0 1 
16 1231 0.01 4.93 0 1 
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TABLE B-2 


DROP-SIZE DISTRIBUTION, SETTLING VELOCITIES Yj, REFLECTION 
COEFFICIENTS Yj; AND IMPACTION EFFICIENCIES Ej FOR 
TRIAL 3 OF THE WITHLACOOCHEE SPRAY TRIALS, 
FLORIDA, 1980 


Drop-Size Mean Drop 
Category Diameter 
j (um) 
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0. 
OF 
Q. 
Os 
0. 
Q. 
Of 
Oe 
Q. 
0. 
0. 
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TABLE B-3 


DROP-SIZE DISTRIBUTION, SETTLING VELOCITIES V;, REFLECTION 
COEFFICIENTS Yj; AND IMPACTION EFFICIENCIES Ej FOR 
TRIAL 5 OF THE WITHLACOOCHEE SPRAY TRIALS 

FLORIDA, 1980 


Drop-Size Mean Drop . \. 
ae Tae F (m 57) Yj 
1 Q. 0. Bo 
2 0. oF 0 
3 Q. 0. 0 
4 0. OF 0 
5 oF. o2 0 
6 0. Le 0 
7 118%. 1m 0 
8 ab. i€ 0 
9 07 id 0 
10 a). de 0 
11 0. ae 0 | 
12 QO. 2 0 
13 Q. 3. 0 
14 0. 3. 0 
15 03 4. 0 
16 0. 4. 0 
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TABLE B-4 


DROP-SIZE DISTRIBUTION, SETTLING VELOCITIES ¥ 
COEFFICIENTS Yj AND IMPACTION EFFICIENCI 
TRIAL 6 OF THE WITHLACOOCHEE SPRAY T 

FLORIDA, 1980 















Drop-Size 
Category 
ius 


Mean Drop 
Diameter 


(um) 











1 53.6 0. 0 
2 110 0.009 0.354 
3 172 0.02 0.686 
4 228 0.03 0.921 
5 279 0.04 1.19 

6 349 0.10 153 

7 427 0.10 1.84 

8 491 0.10 2.10 

9 583 0.20 2.53 
10 682 0.10 3rd 
11 769 0.10 ae 

12 862 0.10 3 il 

13 937 0.04 3.99 

14 1019 0.03 4.29 

15 1100 0.02 4.55 

16 1217 0.01 4.92 


_, REFLECTION 
ES £4 FOR 
RIALS 
Yj 5 
0.45 0.227 
0 0.954 
0 1 
0 1 
0 1 
0 1 
0 1 
0 1 
0 1 
“0 1 
0 1 
0 1 
0 1 
0 1 
0 1 
0 1 


TABLE B-5 


DROP-SIZE DISTRIBUTION, SETTLING VELOCITIES V3, REFLECTION 
COEFFICIENTS Yj; AND IMPACTION EFFICIENCIES Ej FOR 
TRIAL 7 OF THE WITHLACOOCHEE SPRAY TRIAL 








FLORIDA, 1980 





Drop-Size Mean Drop 
Category Diameter i ae 
j (um) : : : 
1 64.7 0.001 males) O31 0.483 
2 120 0.009 0.417 0 1 
3 181 0.02 0.724 0 1 
4 230 0.03 0.942 0 1 
2) 275 0.04 aly, 0 1 
6 337 0.10 1.47 0 i 
7 411 0.10 Liki 0 1 
8 477 0.10 22a 0 1 
9 582 0.20 (age 0 1 
10 694 0.10 3:05 0 1 
11 782 0.10 3736 0 1 
12 898 0.10 eRash. 0 ih 
13 1052 0.04 4.33 0 ] 
14 i222 0.03 4.86 0 1 
15 1330 0.02 ae 0 1 
16 1390 0.01 5% on 0 1 
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TABLE B-6 


DROP-SIZE DISTRIBUTION, SETTLING VELOCITIES Vj, REFLECTION 
COEFFICIENTS Yj AND IMPACTION EFFICIENCIES Ey FOR 
TRIAL 10 OF THE WITHLACOOCHEE SPRAY TRIAL 
FLORIDA, 1980 





Drop-Size Mean Drop 

Category Diameter E; 
j (um) 
1 59.1 0.001 0.149 23 0.115 
2 108 0.009 0.347 0 0.385 
3 161 0.02 0.635 0 0.855 
4 205 0.03 0.824 0 ii 
5 246 0.04 1 01t 0 1 
6 308 0.10 mele35 0 1 
/ 383 0.10 1.66 0 1 
8 444 0.10 1.91 - Q 1 
9 aca 0.20 2.28 0 1 
10 625 0.10 Mees) 5 0 1 
skal n/ 02 0.10 30 0 i; 
12 793 O AG 3.44 0 1 
13 882 0.04 3.78 0 1 
14 997 0.03 4.20 0 ] 
15 1110 0.02 4.56 0 ‘i 
16 1226 0.01 4.93 0 i 
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TABLE B-7 


DROP-SIZE DISTRIBUTION, SETTLING VELOCITIES Vj, REFLECTION 
COEFFICIENTS Y; AND IMPACTION EFFICIENCIES Ey FOR 


TRIAL 11 OF THE _WITHLACOOCHEE SPRAY TRIAL 
FLORIDA, 1980 













Drop-Size Mean Drop 
Category Diameter 
j (um) 









1 51.2 0.001 0.117 46 
2 105 0.009 0.341 0 
3 150 0.02 0.585 0 
4 203 0.03 0.817 0 
5 262 0.04 sa 0 
6 335 0.10 1.46 0 
7 420 0.10 1.81 0 
3 491 0.10 2.10 0 
9 604 0.20 2.64 0 

10 719 0.10 3.14 0 

11 788 0.10 3.40 0 

12 877 0.10 3.73 0 

13 966 0.04 4.06 0 

14 1039 0.03 4.29 0 

15 1151 0.02 4.64 0 

1323 0.01 ar 0 


- 
[@) 
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TABLE B-8 - 


DROP-SIZE DISTRIBUTION, SETTLING VELOCITIES Vj, REFLECTION 
COEFFICIENTS Yj AND IMPACTION EFFICIENCIES E; FOR 
TRIAL 13 OF THE WITHLACOOCHEE SPRAY TRIALS 
FLORIDA, 1980 
















Drop-Size Mean Drop 
Category Diameter 
j (um) 






1 Um 0.001 Q. 

2 LBs 0.009 OF 

3 145 0.02 Q. 

+ 184 0.03 Oe 756 0 
5 220 0.04 0.887 0 
6 ova 0.10 118 0 
7 342 0.10 1.50 0 
8 400 0.10 Lis 0 
9 486 0.20 oOF 0 
10 5/0 OO 2.47 0 
11 637 0... 10. Cd 72 
2 736 0.10 325 0 
13 825 0.04 S256 0 
14 878 0.03 3.26 0 
ale) 928 0.02 3.94 0 
16 1008 0.01 4.23 0 
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APPENDIX C 
DEPOSITION PATTERNS AT THE TOP OF THE SLASH PINE CANOPY 
AND BENEATH THE CANOPY 


Figures C-1 through C-8 show observed and model deposition at the top 
Ometicessasieoimescancpy s.OlulGiais 2, 3, .5,;'6, 7, 10, 11 and 13. of ‘the 
Withlacoochee Spray Trials. Observed and model ground-level deposition 


below the slash pine canopy for these trials are shown in Figures C-9 
through C-16. 
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FIGURE €-1. Observed and mo de } deposition at the top of the 
slash pine canopy for Trial 2, Withlacoochee 
Spray irialse Florida, 1980. 
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FIGURE C-2. Observed and model deposition at the top of the 
slash pine canopy for Trial 3, Withlacoocnee 
Sptayurriars, Florida, 1980, 
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FIGURE C-3. Observed and mode] deposition at the top of the slash pine 
canopy for Trial 5, Withlacoochee Spray Trials, Florida, 
1980. 
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FIGURE C-4. Observed and model deposition at the top of the slash pine 
canopy for Trial 6, Withlacoochee Spray Trials, Florida, 
1980. 
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FIGURE C-5. Observed and model deposition at the top of the slash pine 


canopy for Trial 7, Withlacoochee Spray Trials, Florida, 
1980. : 
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FIGURE C-6. Observed and model deposition at the top of the slash pine 
canopy for Trial 10, Withlacoochee Spray Trials, Florida, 
1980. 
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canopy for Trial ll, Withlacoochee Spray Trials, Florida 
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Observed and model deposition at the top of the slash pine 
1980. 
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FIGURE C-8. 
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Observed and model deposition at the top of the slash 
pine canopy for Trial 13, Withlacoochee Spray Trials 
Por da wei9oul: 
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FIGURE C-9. Observed and model ground-level depositions below the slash 


pine canopy for Trial 2, Withlacoochee Spray Trials, Florida, 
1980. a 
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' FIGURE C-10. Observed and model ground-level deposition below the slash pine 
canopy for Tria] 3, Withlacoochee Spray Trials, Florida, 1980. 
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FIGURE C-ll1. Observed and model ground-level deposition below the slash pine 
canopy for Trial 5, Withlacoochee Spray Trials, Florida, 1980. 
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Observed and model ground-level deposition below the slash pine 
canopy for Trial 6, Withlacoochee Spray Trials, Florida, 1980. 
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Observed and model ground-level deposition below the slash 


pine canopy for Trial 10, Withlacoochee Spray Trials, 


Florida, 1980. 
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FIGURE C-15. Observed and model ground-level deposition below the slash 


pine canopy for Trial 11, Withlacoochee Spray Trials, Florida, 
1980. C-16 
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FIGURE C-16. 
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Observed and model ground-level deposition below the slash pine 
canopy for Trial 13, Withlacoochee Spray Trials, Florida, 1980. 
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